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Abstract 

We investigate a possibility of precision measurements for parameters of the 
Littlest Higgs model with T-parity at the International Linear Collider (ILC). 
The model predicts new gauge bosons (An, Zn, and Wu), among which the 
heavy photon (Ar) is a candidate for dark matter. The masses of these new 
gauge bosons strongly depend on the vacuum expectation value that breaks 
a global symmetry of the model. Through Monte Carlo simulations of the 
processes: e + e~ — > AuZu and e + e~ — > Wj[W H , we show how precisely the 
masses can be determined at the ILC for a representative parameter point of 
the model. We also discuss the determination of the Little Higgs parameters 
and its impact on the future measurement of the thermal abundance of the 
dark matter relics in our universe. 
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1 Introduction 



There is no doubt that the Higgs boson is the most important particle not only for 
the confirmation of the Standard Model (SM) but also for the exploration of physics 
beyond the SM. Quadratically divergent corrections to the Higgs mass term suggest 
that new physics should appear at the scale around 1 TeV. However, electroweak 
precision measurements require that the scale is larger than O(10) TeV in order not 
to conflict with the measurements pQ. This problem is called the little hierarchy 
problem, and many people expect that new physics involves some mechanism to 
solve the problem. 

There are a number of scenarios for new physics involving such a mechanism. The 
most famous one is the supersymmetric scenario. Recently, alternative one called 
the Little Higgs scenario 0, E] has been proposed. In this scenario, the Higgs boson 
is regarded as a pseudo Nambu-Goldstone boson associated with a global symmetry 
at some higher scale. Though the symmetry is not exact, its breaking is specially 
arranged to cancel quadratically divergent corrections to the Higgs mass term at 
1-loop level. This is called the Little Higgs mechanism. As a result, the scale of new 
physics can be as high as 10 TeV without a fine-tuning on the Higgs mass term. Due 
to the symmetry, the scenario necessitates the introduction of new particles such as 
heavy gauge bosons and top partners. 

It is also known that most of Little Higgs models still suffer from severe con- 
straints from electroweak precision measurements due to direct couplings among a 
new heavy gauge boson and SM particles [3]. In order to resolve the problem, a 
symmetry called T-parity is imposed on the models [5]-[7]. Under the parity, new 
particles are assigned to be T-odd (i.e. with a T-parity of —1), while the SM particles 
are T-even. Thanks to the symmetry, dangerous interactions mentioned above are 
prohibited. Furthermore, the lightest T-odd particle is stable and provides a good 
candidate for dark matter. In this article, we focus on the Littlest Higgs model with 
T-parity as a simple and typical example of models implementing both the Little 
Higgs mechanism and T-paritj0 [S]-[E]- Heavy photon plays the role of dark matter 
in this model [8j [9] . 

In order to test the Little Higgs model, precise determinations of properties of 
Little Higgs partners are mandatory, because these particles are directly related to 

x We assume that T-parity is an exact symmetry. It depends on the UV completions whether 
the T-parity is an exact symmetry or not [101 lllj . 
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the cancellation of quadratically divergent corrections to the Higgs mass term. In 
particular, measurements of heavy gauge boson masses are quite important. Since 
heavy gauge bosons acquire mass terms through the breaking of the global symmetry 
mentioned above, precise measurements of their masses allow us to determine the 
most important parameter of the model, namely the vacuum expectation value of 
the breaking. Furthermore, because the heavy photon is a candidate for dark matter, 
the determination of its property gives a great impact not only on particle physics 
but also on astrophysics and cosmology. At the Large Hadron Collider (LHC), 
top partners are expected to be copiously produced, and their properties will be 
determined accurately [T2"| [T3] . However, it is difficult to determine the properties 
of heavy gauge bosons at the LHC, because they have no color charge [Lf] . 

On the other hand, the International Linear Collider (ILC) will provide an ideal 
environment to measure the properties of heavy gauge bosons. The ILC is the future 
electron-positron linear collider for the next generation of the high energy frontier 
physics. At the ILC, electrons and positrons are accelerated by two opposing linear 
accelerators installed in an about 30 km long underground tunnel, and are brought 
into collision with a center of mass energy of 500 GeV-1 TeV. Heavy gauge bosons 
are expected to be produced in a clean environment at the ILC, which enables us 
to determine their properties precisely. In this article, we study the sensitivity of 
the measurements to the Little Higgs parameters at the ILC based on a realistic 
Monte Carlo simulation. In addition, from the simulation results, we estimate the 
capability of the ILC to determine the thermal abundance of the dark matter (heavy 
photon) relics in our universe. 

This article is organized as follows. In the next section, we briefly review the 
Littlest Higgs model with T-parity. Simulation framework such as a representative 
point in the parameter space of the model and the simulation tools used in our study 
are presented in section [31 Details of analysis for heavy gauge boson productions at 
the ILC are discussed in section HJ where we show expected measurement accuracies 
of the heavy gauge boson properties at both ^/s = 500 GeV and 1 TeV. In section 
[3J we will show how powerful the ILC is to determine the Little Higgs parameters 
based on the simulation results. Connection to cosmology from the ILC experiment 
is also discussed in this section. Section [6] is devoted to summary. 
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2 Littlest Higgs model with T-parity 



In this section, we briefly review the Littlest Higgs model with T-parity, in particular 
focusing on gauge-Higgs and lepton sectors of the model. (For general reviews of 
Little Higgs models, see Refs. [TBI fTB| . ) 
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2.1 Gauge-Higgs sector 

The Littlest Higgs model with T-parity is based on a non-linear sigma model de- 
scribing an SU(5)/SO(5) symmetry breaking. The non-linear sigma field S is 

S = e 2in// S , (1) 

where / ~ 0(1) TeV is the vacuum expectation value of the breaking. The Nambu- 
Goldstone (NG) boson matrix n and the direction of the breaking S are 



(2) 



Here, we omit the would-be NG fields in the n matrix. An [SU(2)xU(l)] 2 subgroup 
in the SU(5) global symmetry is gauged, which is broken down to the diagonal sub- 
group identified with the SM gauge group SU(2) L xU(l)y. Due to the presence of the 
gauge interactions and Yukawa interactions introduced in the next subsection, the 
SU(5) global symmetry is not exact, and particles in the n matrix become pseudo 
NG bosons. Fourteen (= 24 — 10) NG bosons are decomposed into representations 
lo © 3o © 2±i/2 © 3±i under the electroweak gauge group. The first two representa- 
tions are real, and become longitudinal components of heavy gauge bosons when the 
[SU(2) xU(l)] 2 is broken down to the SM gauge group. The other scalars 2±i/ 2 and 
3±i are a complex doublet identified with the SM Higgs field (H in Eq. ([2])) and a 
complex triplet Higgs field ($ in Eq. (j2J)), respectively. 
The kinetic term of the S field is given as 



Cv, = ^-Tr 



<9 M £ -iV2{g(WZ+ SW T ) + #'(B£ + SB T )} 



(3) 



where W = WfQ% (B = BjYj) is the SU(2) j (U(l) j ) gauge field and g (g') is the 
SU(2) L (U(l)y) gauge coupling constant. With the Pauli matrix a a , the generator 
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Qj and the hyper-charge Yj are given as 
( a a \ 
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Y x = diag(3, 3, -2, -2, -2)/10 , F 2 = diag(2, 2, 2, -3, -3)/10 . (4) 
It turns out that the Lagrangian in Eq. ([3]) is invariant under T-parity, 

n «-> -nnn, «-> B a «-> s 2 ; n = diag(i, 1, -1, 1, 1). (5) 

This model contains four kinds of gauge fields. The linear combinations W a = 
(Wi + W%)/\/2 and B = (B 1 + B 2 )/V2 correspond to the SM gauge bosons for the 
SU(2)x and U(l)y symmetries. The other linear combinations W£ = (W^ — W£)/y/2 
and B u = {B 1 - B 2 )/V2 are additional gauge bosons called heavy gauge bosons, 
which acquire masses of O(f) through the SU(5)/SO(5) symmetry breaking. After 
the electroweak symmetry breaking with (H) = (0,f/v / 2) T , the neutral components 
of W# and Bn are mixed with each other and form mass eigenstates An and Zh, 



The mixing angle #h is given as 
tan #h 



cos $h — sin #h 
sin #h cos 9u 
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where m n = g 2 f 2 (c 2 + 7)/8, m 12 = gg'f 2 (l - cJ)/8, m 22 = g' 2 f 2 (5c 2 + 3)/40, 
and Cf = cos(\r2v/f). Since the mixing angle is considerably suppressed, A H is 
dominantly composed of B^. Masses of gauge bosons are given by 
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i fmn + m 22 + yj (mn - ™22) 2 + 4mf 2 
i ^mn + m 22 - ^/ (mn - "I22) 2 + 4mf 2 
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(2, -3/10; 1,-1/5) 


1(2) 


(1,-1/5; 2, -3/10) 


eR 


(1,-1/2; 1,-1/2) 



Table 1: Quantum number of [SU(2) x [7(1)] 2 for particles in the lepton sector. 



As expected from the definitions of Ah, Zh, and Wu, the new heavy gauge bosons 
behave as T-odd particles, while SM gauge bosons are T-even. 

Scalar potential terms for H and $ fields are radiatively generated [3l [8], 

V(H, $) = XfTi [$t$] _ ffrfH + j (H^H) 2 + ■■■ . (13) 

Main contributions to fi 2 come from logarithmically divergent corrections at 1-loop 
level and quadratically divergent corrections at 2- loop level. As a result, /i 2 is ex- 
pected to be smaller than f 2 . The triplet Higgs mass term, on the other hand, 
receives quadratically divergent corrections at 1-loop level, and therefore is propor- 
tional to f 2 . The quartic coupling A is determined by the 1-loop effective potential 
from gauge and top sectors. Since both \i and A depend on parameters at the cutoff 
scale A ~ Anf, we treat them as free parameters. The mass of the triplet Higgs 
boson $ is given by m| = A/ 2 = 2m 2 l f 2 /v 2 1 where rrih is the mass of the SM Higgs 
boson. The triplet Higgs boson is T-odd, while the SM Higgs is T-even. 

The gauge-Higgs sector of the model is composed of the kinetic term of the £ 
field in Eq. ([3]) and the potential terms in Eq. (fT31 in addition to appropriate kinetic 
terms of gauge fields W?, Bj and gluon G. It can be seen that the heavy photon 
Ah is considerably lighter than other T-odd particles. Since the stability of Ah is 
guaranteed by the conservation of T-parity, it becomes a good candidate for dark 
matter. 



2.2 Lepton sector 



To implement T-parity, two SU(2) doublets and and one singlet e# are in- 
troduced for each SM lepton. The quantum numbers of these particles under the 
[SU(2)x U(l)] 2 gauge symmetry are shown in Table [TJ With these particles, Yukawa 
interactions invariant under gauge symmetries and T-parity turn out to be 
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where AfW are incomplete SU(5) multiplets, A/" (1) = (/W,0,0) T , AA 2) = (0,0, 



£(») = -a 2 M {n \ and Z (n) = -o 2 {v£\ e[ n) ) T , while tf? is a complete multiplet of 
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5*0(5), i = (l c ,Xi,l c ) T - The indices x,y,z run from 3 to 5 whereas i,j = 1,2. 
For X, there are two possible choices: X = (T.^) -1 ^ 4, and X = (Sg 3 ) 1//4 [17]. With 
£ = SofiS^fiSo and £ = £ 2 £o> these fields transform under T-parity as 

7V {1) «-> -S 7V (2) , «-> -tff, S «-> S, X «-> X, £ ^ fi^fi. (16) 

The linear combination = — l^ 2 ')/y2 gives the left-handed SM lepton, 
which acquires the Dirac mass term with cr in Eq. (fT4l) through the electroweak 
symmetry breaking. On the other hand, another linear combination /h = iy^' + 
l^ 2, )/y/2 acquires the Dirac mass term of O(f) with l c = —a 2 (u c ^ n \e c ^) T in Eq. 
(|T5|) . As expected in Eq. (|T6|) . the heavy lepton / H behaves as a T-odd particle, while 
the SM lepton Ism is T-even. The masses of the heavy leptons depend on Ki, 

m eH = V2mf, m un = ^ — — ^ — Kif ~ y/2nif. (17) 

The lepton sector of the model is composed of the Yukawa interactions above and 
appropriate kinetic terms of above leptons involving gauge interactions associated 
with gauge charges shown in Table [TJ 



3 Simulation framework 

3.1 Representative point in the parameter space 

In order to perform a numerical simulation at the linear collider, we need to choose a 
representative point in the parameter space of the Littlest Higgs model with T-parity. 
Firstly, the model parameters should satisfy the current electroweak precision data. 
In addition, the cosmological observation of dark matter relics also gives important 
information. Thus, we consider not only the electroweak precision measurements 
but also the WMAP observation [T5] to choose a point in the parameter space. 
We have calculated the x 2 -function for observables: 

'<*)(») _ /nW X " 

x 2 



^ — \ obs th 



where 0^ s , O®, and AO^ s are an observed value, its theoretical prediction, and the 
error of the observation for observable 0^ l \ For the observed values, the following 
four observables are considered: the W boson mass (m\y = 80.412 ±0.042 GeV), the 
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/ (GeV) 


m h (GeV) 


A 2 




580 


134 


1.5 


0.5 


m All (GeV) 


mw H (GeV) 


mz 3 (GeV) 


m$ (GeV) 


81.9 


368 


369 


440 



Table 2: Representative point used in our simulation study. 



weak mixing angle (sin 9^ = 0.23153 ± 0.00016), leptonic width of the Z boson 
(Ti = 83.985 ± 0.086 MeV) [E], and the relic abundance of dark matter (tt DM h 2 = 
0.119 ± 0.009) [SUj. On the other hand, theoretical predictions for these observables 
depend on th, "ode. pa—; /, xg rn h . the detailed expressions for the 
predictions, see [2TJ [9]). For theoretical predictions, the fine structure constant at 
the Z pole (a _1 (m^) = 128.950), the top quark mass {m t = 172.7 GeV) (22], the 
Z boson mass (jnz = 91.1876 GeV), and the Fermi constant {Gp — 1.16637 x 10 -5 
GeV™ 2 ) [23J have been used as input parameters. For the / parameter, the region 
/ < 570 GeV, which corresponds to rriA H < mw, is unattractive because the pair 
annihilation of A-r into a gauge-boson pair is kinematically forbidden. 

Using the \ 2 function, we have selected a representative point (f,m h , A 2 ) = (580 
GeV, 134 GeV, 1.5). At the representative point, we have obtained oIdm^ 2 of 1.05. 
Notice that no fine-tuning is needed at the sample point to keep the Higgs mass 
on the electroweak scale [211 02] • The masses of the heavy gauge bosons and the 
triplet Higgs boson at the representative point are summarized in Table [2J It can be 
seen that all the heavy gauge bosons are lighter than 500 GeV, which allows us to 
consider their pair production at the ILC. 

Here, we add a comment on the parameter Ki in Eq. ffl5|) . because cross sections 
to produce the heavy gauge bosons depend on the masses of the heavy leptons as well 
as the other parameters mentioned above. Though the parameter ki is not directly 
related to the observables used in the x 2 - ana lysis, it is also constrained by collider 
experiments. Since small Ki means the existence of light en, too small Ki has been 
ruled out by non-observation of new charged particles. On the other hand, large K\ 
is disfavored because it gives a large contribution to four- Fermi operators [21], [T3] . 
Therefore, ki is expected to be 0(1), and we set k\ = 0.5 in this article^. 



2 Unlike the masses of heavy gauge bosons, those of top partners depend on not only / but also 
A2. Since the top partners are irrelevant to our analysis, we will not discuss the parameter. See 
Ref.[8] for more details. 

3 When K[ < 0.45, heavy leptons en and v^i are lighter than heavy gauge bosons Wh and Z^. 
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e+e~ -> A n Z n 


e + e -> Zh^h 


e+e" -> W+WjT 


500 GeV 


1.91 (fb) 






1 TeV 


7.42 (fb) 


110 (fb) 


277 (fb) 



Table 3: Cross sections for the production of heavy gauge bosons. 




Figure 1: Diagrams for signal processes; e + e — ► A^Zr and e + e — ► 17 + W H . 

There are four processes whose final states consist of two heavy gauge bosons: 
e + e~ — > AftAft, A^Z^, Z-rZ-u, and W + W#. The first process is undetectable, thus 
not considered in this article!. The cross sections of the other processes are shown 
in Table [31 It can be seen that the largest cross section is expected for the fourth 
process, which is open at yfs > 1 TeV. On the other hand, because + m^ H is 
less than 500 GeV, the second process is important already at the y/s = 500 GeV. 
We, hence, concentrate on e + e~ — > A^Zu at y/s = 500 GeV and e + e~ — > W^W^ at 
y/s = 1 TeV. Feynman diagrams for the signal processes are shown in Fig. [lj Note 
that Zh decays into A^h, and W# decays into AhW /± with almost 100% branching 
fractions. 



3.2 Simulation tools 



We have used MadGraph [25J to generate signal events of the Little Higgs model, 
while Standard Model events have been generated by Physsim [26^ . We have ig- 



nored the finite crossing angle between the electron and positron beams. In both 
event generations, helicity amplitudes were calculated using the HELAS library [27J, 
which allows us to deal with the effect of gauge boson polarizations properly. Phase 
space integration and the generation of parton 4-momenta have been performed by 
BASES/SPRING [28]. Parton showering and hadronization have been carried out 



Collider signals will be changed drastically in that case [14] , 

4 Furthermore, even if we consider the process e + e~ — * Ah^IhT, its cross section is strongly 

suppressed due to the small coupling between An and leptons. 

5 Initial-state radiation and beamstrahlung have not been included in the event generations. 
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Detector Performance Coverage 

Vertex detector 5 b < 5 © 10/p/3sin 3/2 9 (//m) | cos 6»| < 0.93 

Central drift chamber Sp t /p^ < 5 x 10~ 5 (GeV/c) -1 | cos9\ < 0.98 

EM calorimeter cr E /E = 17%/ y/E © 1% | cos#| < 0.99 

Hadron calorimeter a E /E = 45%/ y/E © 2% | cos6»| < 0.99 

Table 4: Detector parameters used in our simulation study. 



by using PYTHIA6.4 [29], where final-state tau leptons are decayed by TAUOLA 
[30] in order to handle their polarizations correctly. 

The generated Monte Carlo events have been passed to a detector simulator 
called JSFQuickSimulator, which implements the GLD geometry and other detector- 
performance related parameters [31] . In the detector simulator, hits by charged par- 
ticles at the vertex detector and track parameters at the central tracker are smeared 
according to their position resolutions, taking into account correlations due to off- 
diagonal elements in the error matrix. Since calorimeter signals are simulated in 
individual segments, a realistic simulation of cluster overlapping is possible. Track- 
cluster matching is performed for the hit clusters in the calorimeter in order to 
achieve the best energy flow measurements. The resultant detector performance in 
our simulation study is summarized in Table HI 



4 Results from simulation study 

In this section, we present some results from our simulation study for heavy gauge 
boson productions. The simulation has been performed at y/s = 500 GeV for the 
A H Z H production and at y/s = 1 TeV for the W# W£ production with an integrated 
luminosity of 500 fb _1 each. 

4.1 The AnZn production 

The heavy gauge bosons A H and Z H are produced with the cross section of 1.9 fb at 
the center of mass energy of 500 GeV. Since Z H decays into An and the Higgs boson, 
the signature is a single Higgs boson in the final state, mainly 2 jets from h — > bb (with 
a 55% branching ratio). We, therefore, define A^Z^ — > AnA^bb as our signal event. 
For background events, contribution from light quarks was not taken into account 
because such events can be rejected to negligible level after requiring the existence 
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Process 


uross sec. [tbj 


^ of events 


^ of events after all cuts 


A 1-7 A A n 

A h Zh — >■ A H A H OD 


i.Oo 


k 

ozo 


III 


z/z//i — > z/z/66 


o4.U 


1 ( ,000 


o,ooy 


z7/i — > wbb 


0.0 1 


Z, 1 oO 


1,4UD 


tt -> 


496 


248,000 


264 


zTzT — > 


25.5 


12,750 


178 


z/vZ — > 


44.3 


22,150 


167 


7Z — > 766 


1,200 


600,000 


45 



Table 5: Signal and backgrounds processes considered in the A^Zr analysis. 



of two 6-jets, assuming a 6-tagging efficiency of 80% for 6-jets with 15% probability 
to misidentify a c-jet as a fe-jet. This 6-tagging performance was estimated by the 
full simulation assuming a typical ILC detector. Signal and background processes 
considered in this analysis are summarized in Table [51 Figure [2] shows a typical 
AyiZyi event as seen in the detector simulator. 

The clusters in the calorimeters are combined to form a jet if the two clusters 
satisfy < y cut . y^ is defined as 

_ 22^(1 -oosflg) 
y*i ~~ 77.2 ' v ±a / 

vis 

where % is the angle between momenta of two clusters, E^ are their energies, and 
E vis is the total visible energy. All events are forced to have two jets by adjusting y cut . 
We have selected events with the reconstructed Higgs mass in a window of 100-140 
GeV. In order to suppress the vvh — ► wbb background, the transverse momentum 
of the reconstructed Higgs bosons (px) is required to be above 80 GeV. This is 
because the Higgs bosons coming from the WW fusion process, which dominates 
the vvh — > wbb background, have px mostly below W mass. Finally, multiplying 
the efficiency of double 6-tagging (0.8 x 0.8 = 0.64), we are left with 272 signal 
and 5,419 background events as shown in Table [51 which corresponds to a signal 
significance of 3.7 (= 272/V5419) standard deviations. The indication of the new 
physics signal can hence be obtained at \fs = 500 GeV. 

The An and Zh boson masses can be estimated from the edges of the distribu- 
tion of the reconstructed Higgs boson energies. This is because the maximum and 
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Figure 2: A typical event of A-rZ-r in the simulator, 
minimum Higgs boson energies (E max and .Emm) are written in terms of these masses, 

£max = lZ H E* h + Pz H lZ H P* h , 

E min = lZ H E* h - Pz H lZ H P* h , (20) 

where Pz H (lz H ) is the ^(j) factor of the Zr- boson in the laboratory frame, while 
EhiPh) i s ^ ne energy (momentum) of the Higgs boson in the rest frame of the bo- 
son. Note that E* h is given as (Mf H + M\ - MI h )/(2M Zh ). 

The energy distribution of the reconstructed Higgs bosons with remaining back- 
grounds is depicted in Figj3](a). The signal distribution after backgrounds have been 
subtracted is shown in Figj3](b). The endpoints, E max and E min , have been estimated 
by fitting the distribution with a line shape determined by a high statistics signal 
sample. The fit resulted in m^ H and m Zn being 83.2 ± 13.3 GeV and 366.0 ± 16.0 
GeV, respectively, which should be compared to their true values: 81.85 GeV and 
368.2 GeV. 
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Figure 3: (a) Energy distribution of the reconstructed Higgs bosons with remaining back- 
grounds after the mass cut. (b) Energy distribution of the Higgs bosons after subtracting 
the backgrounds. The distribution is fitted by a line shape function determined with a 
high statistics signal sample. 

4.2 The W+W u production 

H^M^ production has large cross section (277 fb) at the ILC with the center of 
mass energy of 1 TeV. Since decays into and W ± with the 100% branching 
ratio, analysis procedure depends on the W decay modes. In this analysis, we have 
used 4-jet final states from hadronic decays of two W bosons, IfgW^ — ► A^Auqqqq. 
Signal and background processes considered in the analysis are summarized in Table 

El 

All events have been reconstructed as 4-jet events by adjusting the cut on y- 



Process 


cross sec. [fb] 


# of events 


# of events after all cuts 


W n W H A R A R qqqq 


120 


60,000 


41,190 


W + W~ — > qqqq 


1307 


653,500 


678 


e + e~W + W~ — > e + e~qqqq 


490 


245,000 


46 


ev e WZ — > ev e qqqq 


24.5 


12,250 


3,797 


Z H Z H -> A H A R qqqq 


18.8 


9,400 


213 


vvW + W~ — > vvqqqq 


7.23 


3,615 


1,597 


ZW + W~ vvqqqq 


5.61 


2,805 


1,533 



Table 6: Signal and background processes considered in the WgRg analysis. 
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values. In order to identify the two W bosons from W# decays, two jet-pairs have 
been selected so as to minimize a x 2 function, 

X 2 = r c M wl - tr M w ) 2 /a 2 Mw + r c M W2 - *M w ) 2 /a 2 Mw , (21) 

where rec Myt/i(2) is the invariant mass of the first (second) 2-jet system paired as a W 
candidate, tT Mw is the true W mass (80.4 GeV), and ayi w is the resolution for the 
W mass (4 GeV). We required x 2 < 26 to obtain well-reconstructed events. Since 
An bosons escape from detection resulting in a missing momentum, the missing 
transverse momentum ( miss |>r) of the signal peaks at around 175 GeV. We have thus 
selected events with miss p T above 84 GeV. The numbers of events after the selection 
cuts are shown in TableO Notice that the ZnZn and e + e~W + W~ events are reduced 
to a negligible level after imposing all the cuts. The number of remaining W + W~ , 
eu e WZ, vvW + W~ , and ZW + W~ background events is much smaller than that of 
the signal. 

As in the case of the AnZn production, the masses of An and Wn bosons can 
be determined from the edges of the W energy distribution. The energy distribution 
of the reconstructed W bosons is depicted in Figfjja). After subtracting the back- 
grounds from FigJlfa), the distribution has been fitted with a line shape determined 
by a high statistics signal sample as shown in FigJU^b). The fitted masses of An and 
Wh bosons are 81.58 ± 0.67 GeV and 368.3 ± 0.6 GeV, respectively, which are to 
be compared to their input values: 81.85 GeV and 368.2 GeV. Figure [5] shows the 
probability contours for the masses of An and We. at 1 TeV together with that of 
An and Zn at 500 GeV. The mass resolution improves dramatically at y/s = 1 TeV, 
compared to that at y/s = 500 GeV 

The production angle of Wn bosons can be calculated with 2-fold ambiguity from 
the momenta of W bosons, assuming back-to-back production of W# and Wr . It 
turned out that the wrong solutions have a similar distribution to that of the correct 
ones. In Fig® we histogram the two solutions for the production angle, whose 
distribution is consistent with W H being spin-1 particle. 

The angular distribution of jets in the helicity-frame of the parent W carries 
information on the polarization of the W , from which we can extract information 
on the decay vertex of the parent particle. Figure [7] shows the angular distribution 
of the reconstructed jets in the helicity-frame of the reconstructed bosons. The 
distribution indicates the dominance of the longitudinal bosons, implying that 
this coupling arises from the electroweak symmetry breaking. 
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Wenergy(GeV) Wenergy(GeV) 

Figure 4: (a) The energy distribution of the reconstructed W bosons with remaining 
backgrounds after the selection cuts, (b) The energy distribution of the W bosons after 
the subtraction of the backgrounds. The distribution is fitted by a line shape function 
determined with a high statistics signal sample. 
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Figure 5: Probability contours corresponding to (a) 1- and 2-a deviations from the best 
fit point in the An and Zh mass plane, and (b) 1-, 3-, and 5-a deviations in the An and 
Wh mass plane. The shaded area in (a) shows the unphysical region of niA H + mz H > 500 
GeV. 
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Figure 6: Distribution of the production angles of Wn bosons calculated from the mo- 
menta of the reconstructed W bosons. Although there are two solutions, they have similar 
distributions. 
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Figure 7: Angular distribution of the reconstructed jets from W in the W helicity-frame. 
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Figure 8: The simulated cross section measurements for e + e~ — > W^W^ with expected 
statistical errors, which are too small to be seen, as a function of the electron-beam polar- 
ization, where the positron-beam polarization was set to 0%. The minus (plus) sign is for 
the left(right)-handed polarization. 

At the ILC, the electron beam, and optionally the positron beam as well, can be 
polarized. Changing the beam polarization, we can, hence, determine the SU{2)i J 
and U{l)y charges of Wh through the measurements of the e + e~ — > W^W^ cross 
sections. We have studied the measurement accuracy of the cross sections for the 
electron-beam polarizations of -80%, 0%, and +80%, with the positron-beam po- 
larization set to 0%, where the minus (plus) sign is for the left(right)-handed po- 
larization. Figure [8] shows the simulated cross section measurements (data points 
with error bars which are too small to be seen) as a function of the electron-beam 
polarization together with a line representing the prediction by theory. Notice that 
the measured cross sections extrapolate to zero for the 100% right-handed electron 
beam, indicating that Wh has S77(2)l charge but no U(1)y charge. 

5 Discussions 

As shown in the previous section, the masses of the heavy gauge bosons ^h, -^h, 
and Wh can be determined very accurately at the ILC experiment. It is important 
to notice that these masses are obtained in a model-independent way, so that it 
is possible to test the Little Higgs model by comparing them with the theoretical 
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Figure 9: The cross sections of (a) e + e~ — ► A^Z^ and (b) e + e~ — ► W^J" processes as 
a function of . One-sigma regions for the measurement accuracies of these cross sections 
are shown as shaded areas. 

predictions. Furthermore, since the masses of the heavy gauge bosons are from 
the vacuum expectation value / as shown in Eq. ffTU]) . ffTTl) . and ffl~2]) . it is also 
possible to accurately determine /, which is the most important parameter of the 
model. The parameter / is determined to be / = 576.0 ±25.0 GeV from the process 
e+e" -> A u Z n at y/l = 500 GeV, while / = 580.0 ± 0.7 GeV from the process 
e + e~ — > at y/s = 1 TeV. Note that the input value of / is 580 GeV in our 

simulation study. 

Another Little Higgs parameter k\ can also be determined from the results ob- 
tained in the previous section, because production cross sections for the heavy gauge 
bosons depend on the masses of heavy leptons en and v^. Figure [9] shows the cross- 
sections for Ah^h production at ■sfs = 500 GeV and W-^W 7 ^ production at yfs = 1 
TeV as a function of Ki. The measurement accuracies for these cross sections turned 
out to be 22.1% at 500 GeV and 0.8% at 1 TeV, which are shown as shaded regions 
in the figure. These cross section measurements constrain k\. Since the input value 
of Ki is 0.5 in our simulation study, these results correspond to the sensitivity to K\ 
of 9.5% at 500 GeV and 0.8% at 1 TeV. Although there are two possibilities for the 
value of Ki at 1 TeV, we can reject K\ of ~ 0.75 by the measurement at 500 GeV. 

Once we obtain the Little Higgs parameters as above, it is possible to establish 
the connection between cosmology and the ILC experiment. Since the Little Higgs 
model has a candidate for WIMP dark matter [HI [9] , the most important physical 
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0.09 0.1 0.11 0.12 

Figure 10: The probability density of VLh 2 at y/s = 500 GeV and 1 TeV obtained from 
results in our simulation study. The measurement accuracies of cosmological observations 
(WMAP and PLANCK) are also shown as shaded regions. 

quantity relevant to the connection is the thermal abundance of dark matter relics. 
It is well known that the abundance is determined by the annihilation cross section 
of dark matter |32j. In the Little Higgs model, the cross section is determined by 
/ and m h in addition to well known gauge couplings [SJ. The Higgs mass m h is 
expected to be measured very accurately at the ILC experiment [33], so that it is 
quite important to measure / accurately to predict the abundance. 

Figure [TO] shows how accurately the relic abundance can be determined at the 
ILC with the center of mass energies of 500 GeV and 1 TeV. The probability density 
of Qh 2 , which is obtained from the results in the previous section, is depicted. As 
shown in the figure, the abundance will be determined with O(10%) accuracy even 
at \fs = 500 GeV, which is comparable to the WMAP observation. At -y/i = 1 TeV, 
the accuracy will improve to 1% level, which stands up to that expected for future 
cosmological observations such as from the PLANCK satellite [34J . The measurement 
accuracies of these cosmological observations are also shown in the figure in order to 
see the connection between the ILC experiment and cosmology. 

Finally, we add comments on how our results compare with those expected from 
the LHC. Since the LHC is a hadron collider, it is not easy to identify heavy gauge 
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boson productions [H]. However, new colored particles such as top partners will be 
copiously produced and their signals will be detected. When the masses of the top 
partners are around 1 TeV, the Little Higgs parameters such as / will be determined 
with (9(10)% accuracy from the signals [13]. It is then possible to determine the 
properties of the dark matter model-dependently, namely with the use of the relation 
between masses of the top partners and those of heavy gauge bosons. On the other 
hand, when the masses of the top partners are much heavier than 1 TeV, accurate 
determination of the Little Higgs parameters seems difficult at the LHC, though no 
realistic simulation in such a case is performed yet. While the LHC is exploring the 
colored sector of the Little Higgs model, the ILC will move on to cover the weak 
sector, thereby bringing us deeper understanding of the Little Higgs model. 

6 Summary 

The Littlest Higgs Model with T-parity is one of the attractive candidates for physics 
beyond the Standard Model for it solves both the little hierarchy and dark matter 
problems simultaneously. One of the important predictions of the model is the 
existence of new heavy gauge bosons, where they acquire mass terms through the 
breaking of global symmetry necessarily imposed on the model. The determination 
of the masses are, hence, quite important to test the model. In this article, we have 
performed Monte Carlo simulations in order to estimate measurement accuracies of 
the masses (and of cross sections for heavy gauge boson productions) at the ILC for 
a representative parameter point of the model. 

At the ILC with the center of mass energy of 500 GeV, it is possible to produce 
Ah and Zh bosons with a signal significance of 3.7-sigma level. Furthermore, by 
observing the energy distribution of the Higgs bosons from the decays, the masses 
of these bosons can be determined with accuracies of 16.2% for m^ H and 4.3% for 
m Zn . 

Once the ILC energy reaches y/s = 1 TeV, the process e + e~ — > W# W# opens. 
Since the cross section of the process is large, the masses of and A H can be 
determined as accurately as 0.8% and 0.2%, respectively. Using the process, it is also 
possible to confirm that the spin of is consistent with one and the polarization 
of W ± from the decay is dominantly longitudinal. Furthermore, we have shown 
that the gauge charges of the boson could be measured using a polarized electron 
beam. 
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We have also investigated how accurately the Little Higgs parameters can be 
determined at the ILC From the results obtained in our simulation study, it turns 
out that the vacuum expectation value / can be determined with accuracies of 4.3% 
at \/s = 500 GeV and 0.1% at 1 TeV. Another Little Higgs parameter which is 
relevant to the lepton sector of the model, could also be estimated from production 
cross sections. At the ILC with 500 GeV and 1 TeV center of mass energies, Ki could 
be obtained within 9.5% and 0.8% accuracies, respectively. 

Finally, we have discussed the connection between the ILC experiment and cos- 
mology, focusing on the thermal abundance of dark matter relics, which is the most 
important physical quantity for the connection. We have found that the abundance 
can be determined with 10% and 1% levels at y/s = 500 GeV and 1 TeV, respectively. 
These accuracies are comparable to those of current and future cosmological obser- 
vations for the cosmic microwave background, implying that the ILC experiment will 
play an essential role to understand the thermal history of our universe. 
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